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We investigate the possibility that ultra high energy cosmic rays (E > 10 19 eV) are 
related to the distribution of matter on large scales. The large scale structure (LSS) 
data stems from the recent IRAS PSCz rcdshift survey. We present preliminary 
predictions drawn from an anisotropic distribution of sources which follows the 
galaxy distribution. 

1 Introduction 



Ultra high encrgy_cosmic rays (UHECR) are particles with kinetic energies 
above ~ 10 18 cV.El The nature of these energetic particles is presently un- 
known. The reason is twofold: i) These particles interact on the top of the 
Earth's atmosphere producing extensive air showers (EAS) that can be ob- 
served from the ground; in this case, the primary particle can not be observed 
directly; and ii) at these ultra high energies (UHE) the fluxes are extremely 
low (less than 1 particle per square kilometer per year). This makes it im- 
practicable to observe the UHE particles directly using balloons, satellites or 
spacecrafts due to their small acceptance. 

The observation methods are indirect and rely on the observation of the 
secondary particles produced in the EAS. The hadronic particles, as well as 
muons and electrons created by the interactions of the primary particle in the 
atmosphere, are detected on the ground. The EAS also produces detectable 
fluorescent light photons due to the excitation of nitrogen molecules in the air 
by the charged secondary particles. Yet, the secondary charged particles that 
travel with velocities higher than the velocity of the light in the air generate 
Cherenkov radiation that can also be detected. Despite the fact that the 
EAS can be observed by the detection of different kinds of secondaries using 
various techniques, the determination of the nature of the primary particle is 
very difficult and model dependent. As the fluxes are low it is necessary to use 
large ground arrays, and/or many fluorescent light and Cherenkov radiation 
detectors. Until the present moment, the ground arrays and the fluorescent 
light detectors have gathered only a handful of events in the UHE range. 
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The number of events has not been enough to tell us whether the sources are 
cxtragalactic or are located in the Galaxy. 

If the primary particle is a 7-ray or a neutrino the arrival direction would 
point back directly to the source. This would also be the case for charged 
particles if the Galactic magnetic field is < 10~ 6 G and the extragalactic 
magnetic fields are < 1CP 9 G in the case of extragalactic sources. The distance 
to the source is also constrained for most kinds of primaries: If the primary 
particle is a nucleus or a proton (antiproton) the distance to the source is 
limited to less than ~ 100 Mpc for particles with arrival energies above ~ 
6 x 10 19 eV (GZK cutoff, see Fig.(l)). Due to interactions with the cosmic 
microwave background (CMB) photons these particles rapidly lose energy. 
Sources of 7-rays must be even closer because of the short mean absorption 
length for the UHE photons traveling in the CMB. 

The mechanism that provides particles with UHE is also not known. The 
ignorance about the sources makes it harder to determine the mechanism at 
work. For the UHE events no source candidate in the vicinity of the region 
to where the arrival direction points back has been found yet. On the other 
hand, some analysis of the showers profile have been favoring protons as the 
primary particle.B As it was mentioned above, the number of UHE events is 
still too small to allow us, based on the statistics, to answer questions about 
their isotropy and composition. 

In this work, we assume that the UHECR primary particles above 10 19 
eV are predominantly extragalactic protons and that the sources are related 
to the distribution of matter on large scales. It means that without specifying 
the sources themselves or the acceleration mechanism, we would expect an 
excess of events coming from regions with mass overdensitics and less events 
coming from regions with mass underdensities. In the § 2 we briefly describe 
the formalism used; the propagation code is described in the § 3; and the 
smoothing procedure of the density field is described in § [I| Our results are 
presented in the § 5. 



2 Formalism 

In this contribution we apply a generalization of the formalism described in 
Waxman, Fisher and Piran.13 We model the population of UHECR sources S 
within a "box" AV centered at (z, O) as drawn from a Poisson distribution 



S(z,Q) 
ST 



s 



prob s (z, 0) 



exp — S(z, Cl) 
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whose mean value is S(z,Q) = s(z)*B 6p(z,£l) AV. Here s(z) denotes the 
average comoving number of UHECR sources at redshift z and 2$ is some bias 
functional of the local galaxy distribution Sp(z,Cl). The generating function 
{g.f.) of such a distribution is 

f s (u;z,(l) = e X p[s(z,£l)(u-l)], (1) 

where u is a dummy variable. 

The detected number N of UHECR produced by a source within AV with 
observed energy larger than E is also modeled by a Poisson distribution 

prob^ (>E) = N(E N f N exp { -N(E, z) } , 

with mean value 



N(E,z)=AT 



n [E mj (E,z)} (1 + z) 
so iirdLiz)' 



where s = s(z = 0), = 4c 2 Hq (2 + z-2\/l + z) for an — 1 Universe; 

*4 and T denote the detector area and observation time, respectively; Ei n j is 
the energy with which a UHECR observed with energy E was produced at 
redshift z; and no is the number of UHE protons emitted by a source per unit 
time and is assumed to be proportional to dN / dEt n j . We have assumed that 
the source differential spectrum is a power law in energy dN/dEi„j cx E in ^ +1 ^ . 
The g.f. of the last probability distribution is given by 

g N (u;z,E) = exp[N(E,z)(u-l)}. (2) 

Hence, it is straightforward to show from equations (Q) and (|^) that the g.f. 
for the probability of observing a total of N events from AV, with an energy 
larger than E, is expressed by 

F(u; z, ft, E) = exp {s(z, Q) [exp (N(E, z){u - 1)) - l] j . 

The overall UHECR distribution coming from a collection of independent 
volume elements AV has the following g.f.: 



F(u; |J AV, E) = J[F (u; z tl Q t: E) , 

i i 

which for a given line of sight (l.o.s.) can be expressed as an integral over z 
F(u; Q, E) = exp ( / S(z, 0) [exp (N(E,z)(u - 1)) - l] dv\ , (3) 
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In defining A(Q) as 



where dV = c \dt/dz\ d\{z){l + z) 1 dz. 



A(Q) = f ^ S{z,Cl)dV, 
Jo 



one can characterize the distribution of UHECR produced by sources along 
the l.o.s. with energy larger than E by the g.f.. 

G(u; n, E) = — — / S(z, n) exp (N(E, z){u - 1)) dV. 
A(fi) Jo 

From equation ([}]) is then straightforward to show that 

F(u;Q,E) = exp|A(fi) G(u;0,^)-1 |, 
which still is a compound Poisson distribution, although G is not Poissonian. 

3 Propagation code 

The propagation equation takes into account energy losses of the UHE protons 
due to: i) Adiabatic expansion of the Universe; ii) e + e~ pair production; 
and Hi) pion production due to interactions with CMB photons. A proton 
observed at present (z = 0) with energy E must have been produced at an 
epoch z with energy Ei n j = Ei n j(E, z). We assume that the influence of the 
magnetic fields on particles with energies E > 10 19 eV is negligible. Figure 
1 shows the decrease of energy as a function of the distance from the source 
for UHE protons. We note that for a proton to be observed with energies 
above ~6x 10 19 eV the source must be within ~ 100 Mpc from the observer 
irrespective to Ei n j. 

4 Smoothed density field 

The galaxy distribution is estimated from the IRAS PSCz redshift surveyB 
We have computed the smoothed density field on a spherical grid up to 
200 hT x Mpc. The Gaussian-smoothed density field at a grid point n is given 
by 



w sm.n 



(4) 



We have divided the sphere in 72 bins of approximately equal area. Radially, 
the bin size increases in proportion to the IRAS PSCz inter-particlc spac- 
ing [n</>(cz)] -1 / 3 . This smoothing scheme is tailored to keep the shot- noise 
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Figure 1. Propagation of UHE protons in the CMB. The lines represent various proton 
injection energies, Ei n j = 10 7 , 10 6 , 10 5 , 10 4 , 10 3 EeV (from top to bottom respectively). 



uncertainties in the density field roughly constant through out the sampled 
volume. A more detailed analysis of the IRAS PSCz density field can be found 
in Branchini et aln 



5 Results 



We have found that the final results are independent of the cosmological 
parameters (f2, A). Thus, we have used f2 = 1 throughout our calculations 
for the sake of simplicity. For the bias functional *8[5(a;)] we have considered 
2$[<5(x)] = l+8(x). Figure ^presents maps of fluctuations in the mean Cosmic 
Ray intensity, 

5cn(E,n)= 4 ^^j -1, (5) 

for E = (6, 10) x 10 19 eV. In the maps we clearly see the regions from where 
an excess and a deficit of UHECR events is expected following the LSS. The 
specific predictions for futuxe experiments as the Auger project and HiRcs 
will be presented elsewhere.B 
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Figure 2. AitofI projection of Sen, for Ei n j m 60 and 100 EeV (7 = 1.0). The heavy contour 
denotes the zero contour. Dark (light) grey contours denote positive (negative) fluctuations 
equally-spaced at 0.20. The long-dashed line represents the Super-Galactic plane. 
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